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ABSTRACT 

We report on low resolution {R k, 3000) spectropolarimetry of the AO supergiant star HD 92207. This 
star is well-known for significant spectral variability. The source was observed on seven different nights 
spanning approximately 3 months in time. With a rotation period of approximately 1 year, our data covers 
approximately a quarter of the star's rotational phase. Variability in the continuum polarization level is 
observed over this period of time. The polarization across the Hq. line on any given night is typically 
different from the degree and position angle of the polarization in the continuum. Interestingly, H/3 is not 
in emission and does not show polarimetric variability. We associate the changes at Hq as arising in the 
wind, which is in accord with the observed changes in the profile shape and equivalent width of Hq, along 
with the polarimetric variability. For the continuum polarization, we explore a spiral shaped wind density 
enhancement in the equatorial plane of the star, in keeping with the suggestion of Kaufer et al. (1997). 
Variable polarization signatures across Hq are too complex to be explained by this simple model and will 
require a more intensive polarimetric follow-up study to interpret properly. 

Subject headings: polarization - stars: early type ~ stars: emission line, Be - stars: mass-loss - stars: winds, outflows 



1. INTRODUCTION 

Spectropolarimetry is a valuable tool for studying a 
number of important properties of stars. Circular po- 
larization is one of the few direct means of measuring 
stellar magnetic fields (e.g., Babcock 1958). Linear 
polarization, such as from Thomson scattering in cir- 
cumstellar envelopes of hot stars, has been important 
for ascertaining deviations from sphericity in circum- 
stellar media (e.g., the Be stars [Poeckert & Marlbor- 
ough 1976; McLean & Brown 1978; Wood, Bjorkman, 
& Bjorkman 1997] and supernovae [Wang et al. 1996; 
Leonard et al. 2005]). A challenge for this area is that 
net polarizations tend to be small so that large tele- 
scopes are required to obtain sufficiently high quality 



data for sources that are only moderately faint. For- 
tunately, astronomers now have access to several 8- 
10 m class or larger telescopes that are outfitted with 
polarimeters. With such instrumentation, studies of 
polarizations across spectral lines - combining good 
spectral resolution with large telescope collecting ar- 
eas - have become increasingly popular and impor- 
tant for understanding structure in circumstellar media. 
Here we present an analysis of variable linear polariza- 
tion across the Hq line and in the adjacent continuum 
for the AO supergiant star HD 92207 at seven epochs 
over a period of three months. 

This particular star was chosen for several differ- 
ent reasons. It has been studied spectroscopically by 
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Kaufer et al. (1996) who found substantial line profile 
and line equivalent width changes at Hq . These varia- 
tions appear to be related to the well-known discrete 
absorption components (or "DACs") that are com- 
monly seen in UV lines of early-type stars (e.g., Massa 
et al. 1995; Kaper et al. 1996). HD 92207 is not far 
removed from the Luminous Blue Variable stars and 
O supergiants that show variable continuum polariza- 
tions (Lupie & Nordsieck 1987; Taylor et al. 1991; 
Harries, Howarth, & Evans 2002; Davies, Vink, & 
Oudmaijer 2007). Moreover, the relatively late spec- 
tral type of HD 92207 within the early-type class indi- 
cates that Ha acts more nearly as a scattering line than 
a recombination line (Puis et al. 1998), thus making 
the star a prime candidate for exploring variable polar- 
ization across a line dominated by wind emission. 

A description of the observations is presented in the 
following section. An analysis of the data is provided 
in section 3. A discussion of the implications of the 
results is given in section 4. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations of the supergiant HD 92207 were 
obtained in service mode from 2007 January to 2007 
March at the European Southern Observatory with 
FORS 1 (FOcal Reducer low dispersion Spectro- 
graph), mounted on the 8-m Kueyen telescope of 
the VLT. This multi-mode instrument is equipped 
with polarization analyzing optics comprising super- 
achromatic half- wave and quarter-wave phase retarder 
plates, and a Wollaston prism with a beam divergence 
of 22" in standard resolution mode. To perform lin- 
ear polarization measurements, a Wollaston prism and 
a half-wave retarder waveplate were used. The wave- 
plate was rotated in 22.5° steps between 0° and 157.5°, 
taking two 20 s sub-exposures at each of the eight po- 
sitions. 

The GRISM600R was used in the FIMS (FORS 
Instrumental Mask Simulator) observing mode in the 
wavelength range 4672-6795 A to cover the hydrogen 
Balmer lines H^ and H^. The spectral resolution of the 
FORS 1 spectra taken with this setting and a 0.4" slit 
was R~3000, corresponding to a velocity resolution of 
about 100 km s"^ 

The readout time was reduced to about 40 s by win- 
dowing the CCD, and the use of a non-standard read- 
out mode (A, 1x1, low) provided a broader dynamic 
range, allowing us to increase the S/N of the indi- 
vidual spectra. Due to the brightness of the target 
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Fig. 1 . — Ha spectra plotted as normalized to the local con- 
tinuum and shifted in the vertical direction for clarity. The 
epochs are indicated for each line. Note that the data have 
been smoothed. Substantial variations in profile shape are 
evident, with line equivalent widths given in Tab. [T] The 
peak emission clearly shifts between epochs. 



(my = 5.5) and excellent seeing conditions, a few 
exposures obtained during the second and the last ob- 
serving night were saturated. For these phases, only 
measurements at the retarder waveplate positions be- 
tween 0° and 67.5° are available. 

For all seven nights we calculated the values 
Pq = Q/I and Pu = U/I, where /, Q, and 
U are the Stokes parameters as defined in Shurcliff 
(1962). The total linear polarization is given by 

P = 



' Pq + P^, and the polarization position an- 
gle is ip = 0.5 ta,n^^ (Pu/Pq). The associated errors 
for the determination of Pq, Pjj, P, and ip are esti- 
mated from error propagation, based on pure photon 
noise in the raw data (see Tab.[T]i. 

The spectropolarimetric calibration of FORS 1 
was checked with the spectropolarimetric standard 
NGC 2024-1. On 2007 January 8, 300 s exposures 
were taken at retarder waveplate angles 0°, 22.5°, 
45°, and 67.5°. These data were reduced in the same 
manner as those obtained for our target HD 92207. 
We obtained the values P = 9.65±0.12% and ip = 
136.82±0.34°, which are in very good agreement 
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Table 1 : Observing Journal 



Date 


MJD 


Ha EW 


Qcont 


Tj-a 
^ cont 






(A) 


(%) 


(%) 


2007 Jan 06 


54106.356 


-A.02 


-1.96 ±0.04 


-2.55 ±0.04 


2007 Feb 24 


54155.189 


-3.16 


-2.12 ±0.04 


-2.60 ±0.04 


2007 Feb 26 


54157.112 


-2.91 


-2.14±0.04 


-2.64 ±0.04 


2007 Feb 27 


54158.117 


-2.60 


-2.13 ±0.04 


-2.63 ±0.04 


2007 Mar 1 1 


54170.291 


-3.50 


-1.90 ±0.04 


-2.51 ±0.04 


2007 Mar 25 


54184.162 


-3.47 


-1.79 ±0.04 


-2.60 ±0.04 


2007 Mar 26 


54185.151 


-2.96 


-1.87 ±0.04 


-2.55 ±0.04 



" The continuum polarizations in the final two columns are taken from the continuum intervals of 6460 — 6470A and 

6630 — 6650A adjacent to the line. 



with the previous measurements P = 9.53±0.02% 
and V = 136.75±0.16° by Fossati et al. (2007), who 
used numerous spectropolarimetric observations of 
NGC 2024-1 retrieved from the ESO archive. 

3. ANALYSIS 

Here we describe the characteristics of the ob- 
served polarizations and their implications for the 
source. Based on previous studies of H^, it is not 
surprising that even in low resolution spectra, the P 
Cygni emission line of of HD 92207 shows sig- 
nificant variability, both in line shape and equivalent 
width. Figure [T] shows average line profile shapes for 
the seven spectra with dates labeled. The different 
spectra have been shifted vertically for better display. 
The equivalent widths (EWs) were evaluated between 
±275 km s^^, the terminal speed Voo of the wind (see 
Tab. |2]l, and are listed in Table [T] Note that all of the 
lines show a net emission. 

Figure |2] shows the polarization of and the 
neighboring continuum in the form of Q-U diagrams 
for each night. Each panel displays the date, mov- 
ing chronologically clockwise from lower left. The 
black point in each panel signifies the mean contin- 
uum polarization from relatively "clean" regions of 
the spectra that appear to lack lines, including 6545- 
6557 A on the shortward side of and 6569-6580 A 
on the longward side. We determined a typical stan- 
dard deviation in the continuum polarization to be ap- 
proximately a « 0.04% for any given night, which 
is approximately the size of the plotted black squares. 
Combining all seven nights, the average polarizations 
are (Pg) = -1.99% and {Pu) = -2.58%; however, 
no attempt has been made to correct the data for inter- 



stellar polarization, so these values cannot be consid- 
ered intrinsic to the star. However, our analysis focuses 
on variable polarization which is intrinsic to the star 

For wavelengths within the line, dashed and dotted 
line types distinguish between blue and red shifts from 
line center (these appear as blue and red colors in the 
online version of the figure). Again, points considered 
to be within the line are those between ±275 km s^^ 
of line center. 

The central panel is an overplot of the line data 
for the other seven panels, indicating the full range of 
polarimetric variations across Hq,. There is a change 
in the continuum polarization over the 3 month time 
span. The overall amplitude of change is about 0.5%, 
which is over 10 times larger than the dispersion for 
any one of the averages. 

Most remarkable are the observed variations in the 
polarization across the Hq P Cygni line. The average 
of the polarization within the line is typically greater 
than the dispersion in the continuum knot by a factor 
of 4 or more. The two main exceptions are the nights 
of March 1 1 that shows very little variation across the 
line and March 26 that shows the greatest variations, 
essentially a change in polarization of about 2%. 

It is notable that the variations across H^ imply dif- 
ferent polarization position angles on different nights. 
Recall that tan {2ip) — Pjj/Pq- Consequently, for 
a fixed geometry, ip will be a constant, implying that 
Pu oc Pq which amounts to variations along a line 
in Q-U diagrams. Although this seems approximately 
true for the wavelength-dependent polarization across 
Hq on any given night (February 26 being an ex- 
ception), the extensions of the polarization from the 
knot of continuum points appear to fall along differ- 
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Fig. 2. — The variable polarization across the line is displayed here in separate panels with dates given. The center panel 
overplots the various QU-loops. The line types are dashed for wavelengths blueward of line center, and dotted for those redward. 
(The online version shows blue and red colors corresponding to the velocity shift from line center.) The black points are the nightly 
mean Q-U values of the nearby continuum. Notably, variations between the polarization in the line and in the continuum is generally 
present, but to varying degrees, and also at different position angles. 



ent orientations between nights. Take for example the 
consecutive nights of March 25 and 26. Polarimetric 
changes are in U for the former, but in Q for the latter, 
suggesting a rotation in ip by 90° after just one day. 
Given that the size and rotation of the star suggest a 
rotation period of approximately 1 year (see Tab. |2]l, 
variations at the level of 1 day seems far too short to 
arise from a substantial change in geometry. 

4. DISCUSSION 

We have presented relatively medium resolution 
spectropolarimetric data of the highly variable A su- 
pergiant HD 92207. It has been suggested by Kaufer 
et al. (1997) that the observed photometric and Ha Hne 
variations are the result of a corotating structure in the 
wind, which they consider to be in the star's equato- 



rial plane. To explore that possibility, we begin with a 
consideration of the variable continuum polarization. 

We have developed a phenomenological model for 
Thomson scattering from an enhanced scattering re- 
gion in a spherical wind in the form of a spiral pattern, 
similar in spirit to the work of Brown et al. (1994) who 
sought to model DACs in emission lines. Our objec- 
tive is to obtain a reference model for the global wind 
morphology that we can use for interpreting the Hq, 
polarization. 

We consider a corotating stucture as a simple spi- 
ral that is top-bottom symmetric about the plane of the 
star's rotational equator Hence, the spiral structure 
has a guiding center that is always in the equatorial 
plane. This center obeys an equation of motion for 
the radial wind velocity law and conservation of an- 
gular momentum in the rotating frame, so it follows a 
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Table 2: Model Parameters 



Factor 


Value 




223i?0 


Vrot sini° 


30kms^i 




1.3 X 10-6 Mq yr-i 




275kms-i 




1 


b" 


0.98 


V' 


25 




55° 


ro 


8° 


*8 


70° 



' Przybilla et al. (2006); Kudritzki et al. (1996); " best fit 
model parameters (see text). 



"streak line" (e.g., see Ignace, Bjorkman, Cassinelli 
1998). We adopt a standard "/3-law" for the radial 
flow: 



(1 - buf , 



(1) 



where u — Ri,/r and 6 < 1 determines the radial 
speed of the wind at its base, with vq — Vqo (1 ~ ■ 

Velocity laws with (3=1 and (3 = 2 were consid- 
ered, and the case of /5 = 1 produced a better match to 
the data. In this case the location of the guiding center 
is given analytically with azimuth by 



Woo 

w{u) 



bin 



w{u) 



(2) 



where wq = vq/voc and 51 = 2-k/ P. 

The density for the spiral-shaped perturbed region 
is treated as an excess of density above the otherwise 
spherically symmetric wind. This density excess in 
the spiral pattern is taken to scale with the spheri- 
cal wind density in our "toy" model, thus we con- 
veniently parametrize the excess by a constant factor 
r] = noxccss/?^sph, for Us-ph the spherical wind den- 
sity. In addition to the solution for the guiding cen- 
ter, we also need the cross-section of the spiral. The 
cross-section (i.e., the intersection of the spiral pat- 
tern with a spherical shell) is treated as circular This 
spherical "cap" is axisymmetric, and so assuming the 
electron scattering is optically thin, the polarization 



from any given slice of the spiral is given by Brown 
& McLean (1978), along with the finite star depolar- 
ization correction factor of Cassinelli, Nordsieck, & 
Murison (1987). Summing up contributions from all 
the caps yields the polarization from the structure as 
a function of rotational phase and viewing inclination 
for the rotation axis of the star. 

Note that our model accounts for occultation of 
scattered light by the intervening star, but only in an 
approximate way. We consider a slice as entirely oc- 
culted if its guiding center lies behind the star, and un- 
occulted otherwise. 

The principal model parameters are the density ex- 
cess rj, the half-opening angle of the spiral 5, an orien- 
tation angle between the observer (5— ?7 axes and those 
of the star system -00, and finally the inclination of the 
rotation axis of the star zo- Using a reduced chi-square 
evaluation for a grid of model polarization light curves, 
Table|2]lists the model parameters that provide the best 
fit to the observed continuum data. The star and wind 
properties of HD 92207 are taken from Przybilla et al. 
(2006), except that the mass-loss rate is taken from Ku- 
dritzki et al. (1999) and does not account for clumping. 
The most reasonable match to the observed continuum 
polarizations in the neighborhood of is shown in 
Figure [3] The upper panel is for Pq and the lower 
one for Pjj, displayed as percentage polarizations. The 
rotational phases depend on the star's rotation period. 
Given the radius and minimum rotation speed from Ta- 
ble|2l the maximum period is Pmax ~ 376 d. The true 
period is Prot = Pmax/ sinzQ, where io is constrained 
from our model fitting. As the ephemeris is not known, 
we assigned the rotational phase "0" to the date of our 
first observation, and the phases appearing in Figure [5] 
represent values for our best fit model at io = 70°. 

For the model fitting, there are seven free parame- 
ters: the four listed above plus three offsets - one for 
rotational phase, a vertical offset in Q, and an indepen- 
dent one in U . With fourteen data points in total, the 
reduced chi-square for our best simultaneous fit to the 
Pq and Pu lightcurves is 1.6. That value is primarily 
a reflection of one discrepant data point, the first one 
in the observed Pu ligthcurve. Our model is inher- 
ently smooth, whereas the line variability indicates the 
presence of variable wind structure. For the continuum 
polarization, this wind structure can act as a source of 
"noise". Better temporal and phase coverage is needed 
to interpret the polarimetric lightcurve more fully. 

The non-zero variation of Py indicates that iq is not 
edge-on to the system, since no APij could result for 
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Fig. 3. — Model Pq (top) and Pu (bottom) light curves 
plotted for two rotational cycles based on the model de- 
scribed in the text with parameters listed in Tab.|2l Observed 
data are plotted as points with error bars. The data points 
have been shifted horizontally and vertically to provide the 
best fits. 

such a perspective when the circumstellar scattering 
environment is top-bottom symmetric. However, APjj 
is much less than APq, suggestive that iq is not close 
to pole-on. Since the interstellar polarization is not 
known, there is freedom to shift the observed points 
with respect to the model in both phase (i.e., to) and 
vertical offset to best match a model curve. The ver- 
tical offsets for Pq and Pjj are independent, but the 
phase shift must be the same for both. There is also 
freedom to rotate the model system relative to the ob- 
server Q-U system, but we found the angle to be small, 
with tpQ = 5°. 

The match between the model and the continuum 
data appears to be best around an inclination of io = 
70° with a half-opening angle of 6 = 55°. A 55° 
opening angle coiTesponds to the Van Vleck angle and 
yields the maximum polarization when other param- 
eters are held fixed, thus allowing us to minimize the 
value of the excess density. Because the spherical wind 
is relatively thin in electron scattering (tc « 0.014 for 
solar abundances with ionized H and He neutral), a 
fairly large value of 77 = 25 is needed to match the 
observed APq and APy. Such a large excess would 



need to be justified through physical modeling of the 
stellar wind and atmosphere. 

Overall, the match is reasonably good except for 
the first data point in the lower panel for Pjj which we 
cannot reconcile to the model after exploring a range 
of opening angles, viewing inclinations, and rj values. 
Based on the work of Davies et al. (2005) for lumi- 
nous blue variables, wind inhomogeneities are likely 
contributing a random contribution of unknown am- 
plitude to the observed polarizations from HD 92207, 
which makes the interpretation of any globally coher- 
ent variation with phase challenging for such sparse 
sampling. Even though our spiral model is simplis- 
tic, it does capture the essential points inherent to any 
corotating structure in terms of the generic trends of 
Pq and Pij with phase, viewing inclination, and den- 
sity. It is certainly clear that much better sampling of 
the polarimetric light curve in the continuum is needed 
to critically assess the existence of a corotating region 
in the wind. 

Turning to the Hq, line, the P Cygni morphology 
indicates that the line is dominated by the stellar 
wind. In contrast, the Hp line is in absorption and 
shows no variable polarization. At the spectral type 
of HD 97702, the n = 2 level of hydrogen can act as 
the effective ground state (Puis et al. 1998) explaining 
why it shows such a strong P Cygni shape. It also 
means that line scattering polarization (e.g., Hamilton 
1947) may influence the observed variable polariza- 
tion in Hq. For example, such effects were explored 
by Jeffery (1987, 1989) in the case of SN1987A. 

Figure m displays the polarized Pq and Pu spectra 
in velocity along with the P Cygni profiles to empha- 
size where polarization changes occur within the line 
in relation to the absorption trough and emission peak. 
The vertical scale for the measured percent polariza- 
tions are shown on the right axes, whereas the scale for 
the continuum normalized P Cygni lines is displayed 
on the left axes. For the first 5 panels, it is typically 
the case that either Pq, Pij, or both vary most around 
the line core over an interval that is about a third of 
the wind terminal speed. One clear exception is for 
March 26 that shows the strongest change in polar- 
ization across the line in Q and aligns quite well with 
the redshifted emission peak, although it appears that 
Pjj actually increases in absolute value. The observed 
polarizations are negative, so zero polarization is to- 
ward the top of each panel. Curiously, there is a depo- 
larization at the same location the night prior, except 
in U instead of Q. The overall depolarizations across 
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Fig. 4. — An overplot of the P Cygni lines against the Pq (dotted curve, or red for the online version) and Pu (dashed curve, 
or blue for the online version) spectra using the format of Fig. |2]to highlight the velocity locations of features. The line is 
normalized to the local continuum, and plotted as the solid curve using the scale on leftside axes; polarizations are shown as percent 
using the scale on the rightside axes. 



the emission peak in the line relative to the continuum 
would seem to be a classic "line effect". Polarization 
is normalized at each wavelength to the total emission. 
If the line formation leads to unpolarized radiation that 
is largely unscattered by free electrons before emerg- 
ing from the wind, then the additional line emission 
tends to reduce the polarization at those wavelengths 
relative to values outside the line. In fact, this can be 
an excellent way of placing an upper limit to the inter- 
stellar polarization to the star. However, the line 
of HD 92207 is perplexing because its behavior is not 
consistent. Moreover, night-to-night variations in the 
Hq, polarization are hard to understand in terms of the 
spiral structure that we have considered for the contin- 
uum polarization given that the star's rotation period is 
of order a year. 

It is worth pointing out that Harries (2000) mod- 
eled wavelength-dependent polarization across H^ for 



the O supergiant C, Pup. His models did not include 
contributions to the polarization from line scattering; 
rather, he obtained changes of polarization across the 
line owing to the influence of line opacity for the po- 
larization produced by Thomson scattering. 

The application of Harries' models to HD 92207 is 
unclear There is the likelihood of a significant line 
scattering contribution to Hq because of HD 92207 's 
much later spectral type as compared to C Pup- Other 
effects such as dynamic changes in the wind or vari- 
able steflar iflumination (c.f., Al-Malki et al. 1999) 
will be needed to interpret the polarization across H^. 
For example, it would be useful to combine the meth- 
ods of Li et al. 2000 and Davies et al. 2007 for po- 
larimetric variability from wind clumping and electron 
scattering with the resonance scattering polarization 
effects explored by Ignace, Nordsieck, & Cassinelli 
(2004) for stationary winds. With limited coverage 
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of the rotational phase and only modest spectral res- 
olution, the existing dataset is too poor to undertake a 
detailed calculation of the Hq, line to interpret the ob- 
served polarizations. With future data, such an effort 
would be worthwhile because the continuum polariza- 
tion constrains the global wind morphology whereas 
the H„ line is sensitive to vector velocity flow. 
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